function [ M ] = fuckers( p, q )
%FUCKERS Summary of this function goes here
%   Detailed explanation goes here

r = p - q;

rn = norm(r);

M = [-3*r(2)*(rn^2+r(1)^2), r(1)*(r(1)^2+r(3)^2-2*r(2)^2), -3*r(1)*r(2)*r(3);
     r(1)*(r(1)^2+r(3)^2-2*r(2)^2), -r(2)*(rn^2+3*r(2)^2), r(3)*(r(1)^2+r(3)^2-2*r(2)^2);
     -3*r(1)*r(2)*r(3), r(3)*(r(1)^2+r(3)^2-2*r(2)^2), -3*r(2)*(rn^2+r(3)^2)] / (16 * pi * rn^5);

end

